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Abstract 

After carefully regularizing the Wheeler - De Witt operator, which 
is the Hamiltonian operator of canonical quantum gravity, we find a 
class of exact solutions of the Wheeler - De Witt equation. 
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The enigma of quantum gravity is probably the most challenging prob- 
lem of modern theoretical physics (for the recent reviews see [jl]]) It is for the 
difficulty of the problem that in spite of its importance, very little progress 
has been made so far. At the moment there are two major approaches to 
quantum gravity. On the one hand, it is being argued that gravity is a low 
energy approximation of some drastically different fundamental theory; the 
major examples of this kind of ideas are supergravities and, more recently, su- 
perstring theories. If one accepts such a proposal, "quantization" of classical 
general relativity does not make much sense. 

There is however another approach in attempt to quantize the classical 
Einstein general theory of relativity. Nowadays the main line of attack lies in 
using canonical quantization technique. In recent years a great deal of excite- 
ment was risen by introduction of Ashtekar variables which were expected 
to simplify the technical problems of older canonical techniques. However, it 
seems that these hopes have not been fulfilled. Even though the use of the 
loop variables solves part of constraints automatically, it turned out that the 
remaining ones are as difficult to solve as the constraints of the traditional 
canonical approach, based on the Wheeler - De Witt equation (given the 
result of the present paper rather more difficult.) Moreover, due to the fact 
that existing solutions are expressed in terms of very exotic variables, their 
physical meaning is rather obscure. This is the reason why we choose to 
analyse again the original approach to quantum gravity which dates back to 
seminal papers of Wheeler and De Witt [@], It should be stressed how- 
ever, that it seems likely that all these three approaches may reveal three 
different regimes of the same ultimate theory (see |jl|] for extensive discussion 
of this point). 

The Wheeler - De Witt approach to quantization of general relativity is 
based on the following set of equations describing "the wavefunction of the 
universe" : 

V^p"^^(x)^(/i„„) =0 (1) 
HvKDH/^ = (-^^P''P"" + (2A + i?)) ^ = 0, (2) 
where n = (IGnG)^^, G is the Newton's constant. 
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p are momentum operators related to the three metric hmn, 

= z6f6l6'{x-y), (4) 

R is the three dimensional curvature scalar, and A the cosmo logical constant. 
We assume also that the three-space is a compact smooth manifold without 
boundaries. 

Equation (||) can be easily solved by assuming that the wavefunction is 
a function of integrals over space of scalar densities (local functionals). We 
will not use here any minisuperspace approximation so we take the lowest 
order local expressions that are diffeomorphism invariant. In what follows 
we will assume that 

^ = ^(V,7^), (5) 

where 

V = J /h(x)d^x, ^ = / \^h{x)R{x)d^x. (6) 

However, as it stands, equation (0) is meaningless. First the functional 
derivatives act at the same point, so acting on any local functional, the first 
term in (Q) produces 5(0). This problem is resolved by regularization - in 
this paper we use the heat kernel to regularize divergent expressions. At 
this point one important remark is in order. The remarkable property of the 
metric representation is that the metric operator acts as multiplication. This 
means that we do not need to introduce any background metric for regular- 
ization purposes (and later work hard to make the final results background 
independent.) What we have for granted are eigenvalues of the metric oper- 
ator on any state. Strictly speaking, in the heat kernel we should use metric 
operators, but these can be identified with their eigenvalues. 

The second problem is related to the choice of ordering of the operator. 
The use of additional freedom given by the ordering parameter a will be 
essential in our procedure. 

Thus we rewrite eq. (0) as 

-«//y&.„„(.)A-(-,!/;*)j^j^- 

J ohki[y) oK^n[x) J 

(7) 
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where K{x, y; t) is a heat kernel satisfying the equation 

(x, t) = y\,)K{x, y- 1) + iR{x)K{x, y; t) (8) 
with the initial condition 

\i^K{x,y;t) = ^'''^^y\ (9) 

Let us note that in (|^) we cannot take a more general ordering involving 
second functional derivatives of since these will again produce (5^^-'(0). 
Equation (H) can be solved pertubatively in powers of t. To the order which 
will be of interest in the present context, the solution reads @] 

^ exp - l^Rmnix)) A-A^ 

K[x, y; t) = 



(47rt) 



1 + ^ - ^) R{x) + 1' Q (e - 1) (10) 



where A™ = — y™. In order to satisfy the Wheeler-De Witt equation we 
have to cancel OR in the above expression so we choose ^ = |- 
We will use below a shorthand notation 

Kkimn{x, y; t) := gkimn{x)K{x, y; t) (11) 

We need one more technical result, namely the value of 

This value cannot be derived by taking derivatives of terms in the expansion 
for K above (it involves terms of all orders in t), but can be found by solv- 
ing perturbatively the equation for ^f^. After straightforward but tedious 
calculations one finds (to the leading order in t) that the expression (^) 
equals 

(rhihmnix) + r^{h2Rmn{x) + hhmn{x)R{x)) + 0{t^)) . 

(47r)2 ^ ^ 
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where 

= ! + = 6, = «!-M_i (13) 

4 3 6 3 12 8 ^ ^ 

Now we are ready to present our main result. We assume that \I/ depends 
only on V and TZ. Substituting this form to the Wheeler - De Witt equation 
(|^ we find terms proportional to 
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d^yKkimn{x,y;t) — — = - — Jh{x)K{x,x;t) (14) 
dhki{x)dhmn{y) 8 ^ 



/ = -'-^) (15) 

,3 ^Kkimn{x,y;t) 6V 
a y 

1 

(4^ 



1_ U-lh, + rHh + 3h)R{x) + Oit-^)) . (16) 
7r)2 V z / 



6hki{x)6hmn{y) 



\fh{x) (^_3^^(^lM + (3^ + 7/8)i?(x)ir(x, X- 1) j (17) 

/ d'yK,i^n{x,y-t) ^f^. .J'^. , = fi?„„(x)i?"^'^(x) - ^i?2(a;)) 

J dhki[x) dhmn[y) ^ \ 8 J 

(18) 

/ci3,K..„(.,,;t)^^ - (-^^(^)) (19) 

^ 6Kkimn{x,y;t) 5TZ _ I /^^_3L 



(5/lA;/(y) Shmn{x) (47?) t V 2 

((62 + h)R\x) - 62i?n.„(x)i?™"(x)) + 0(t^)) . (20) 



where bi are given by ([T3|). 

Now we face the problem as to if renormalize the equation (i.e., replace 
the singular terms involving inverse powers of t) or to find a solution of 



4 



the equation and renormalize the so obtained wave function. Let us make 
this first choice. We use the method proposed by Mansfield which, after 
multiplying by arbitrary function of t and analytical continuation results 
effectively in replacing the powers t~ 2 by pth derivatives of arbitrary function 
0(s) at s = 0, (s = v^) which we denote by These numbers are therefore 
renormalization constants. 

Collecting all terms we have four equations (the coefficients multiplying 
B?{x), RmnR"^"^, R{x) and 1 respectively): 

, |i,<.. f _ (22) 



97^2 dn"^ V 80 

— -abi - - + 7^0^^ 7^77 + 7Ta&3 + 7Ta&2 + /i^^ = 



AdTZdV dn^ V2 ' 8/ V80 2 _ , 

(23) 

' ^ h + 7^77 -^^ h 2A/i2$ = (24) 



8 9V2 V8 2 V m A 

These equations are linear differential equations with constant coefficients 
so the solutions are combinations of exponential functions of V and TZ. 

We will present an explicit solution when the wavefunction depends only 
on V (the general case is not difficult but the equations are rather lengthy). 
In this case only two of the equations survive and we substitute 

^(V) = Ae'^'^, (25) 

Then we get from ( ^31) and (|2^) that uji satisfies 



+ y + 2A.^ - ^J|L (^,...., + . (26) 
The dependence of the solution on /i^ ~ G^^ can be written as 



uJi = -B± JB^ - Cfi^ (27) 



where the constants B and C can be extracted from eq. (pBl). Depending 
on the values of the constants we get real or complex wave functions - the 
interpretation in these cases would be quite different. 
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In the most general case when the function depends on both V and TZ we 
get a relation between the cosmological constant and the Newton's constant 
expressed via the renormalization constants. 



We conclude this paper with a number of comments: 

1. It is natural to ask if one can find a solution depending on higher order 
local functionals like e.g. 5*2 = J{bR'^ + cR„mR™'^)- The procedure 
should be as follows. First, in order to avoid terms of order in the 
heat kernel expansion, we must choose h and c such as to avoid fourth 
order covariant derivatives in the second functional derivative of 5*2 
(otherwise we would have a large number of independent third order 
tensor structures rather impossible to match). This gives one relation 
between the coefficients. Then one should try to match terms with 
independent combinations of curvature and metric tensors up to the 
order three. It is unclear at the moment if this procedure leads to any 
nontrivial solution. 

2. There are several directions of obvious generahzations of our result that 
are currently under investigations. First, one may consider inclusion 
of matter fields. The case of a scalar field and the closely related 
analysis of the cosmological constant problem will be the subject of 
the forthcomming paper. Second, it is of outmost importance for e.g., 
black holes physics to consider a theory based on a compact manifold 
with boundaries. This problem is, however, much more complicated. 

3. It should be stressed that the wave function alone does not provide 
us with much insight into the physics of the solution. We certainly 
lack one important ingredient which is an inner product. In context of 
the solutions presented above, our goal will be to write down the path 
integral measure ii{hmn)dhmn in the form 

dhijvdR j{h^,v,n), 

where Sh^^^ are variations of the metric which leave V and TZ invariant: 
and then integrate over h^j,^ to obtain 
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Only after succeeding this program we will be able to address questions 
concerning hermiticity of operators and their expectation values. 

Alternatively, we may make use of the quantum potential approach 
to quantum gravity which provides us with the four dimensional dy- 
namics related to a given wave function The major virtue of this 
approach is that most physical predictions concerning the evolution do 
not depend on the inner product. 

4. Given the regularized WDW operator we can analyse the issue of 
anomalies. This problem has been investigated in 0, but we found 
ourselves in disagreement with the results of the papers. There are 
also additional, potentially anomalous terms resulting from the general 
ordering chosen in our approach. We are currently investigating this 
problem, and the result will be published soon. 

5. One may wonder why such simple solutions were never found in spite 
of many papers on semiclassical expansion of the WDW equation. The 
reason is that in the semiclassical expansion, as a rule, one employs 
an expansion of the logarithm of the wave function in powers of G~^. 
However, as our result shows, the exact solution is of the form e.g. e^^^ 
with uji being a nonpolynomial function of (eq. (P?])) therefore such 
an expansion is not finite (one should recall that the renormalization 
coefficients 0'-"^ are dimensionful and therefore can mix with powers of 
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